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In this paper, we investigate in general how thermodynamic quantities such as the polarization, 
magnetization and the magneto-electric tensor are affected by the boundaries. We show that when 
the calculation with periodic boundary conditions does not involve a Berry's phase, the quantity in 
question is determined unambiguously by the bulk, even in the presence of gap less surface states. 
When the calculation involves a Berry's phase, the bulk can only determine the quantity up to 
some quantized value, given that (i) there are no gapless surface states, (ii) the surfaces do not 
break the symmetries preserved by the bulk, and (iii) the system is kept at charge neutrality. If 
any of the above conditions is violated, the quantity is then determined entirely by the details at 
the boundary. For components of response functions such as the magneto-electric tensor that are 
entirely determined by the bulk, the value of those components should connect smoothly to finite 
frequency and momentum. 
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I. INTRODUCTION 

Recently there has been a resurging interest in un- 
derstanding the general orbital magneto-electric (ME) 
responsei— . This is due to the fact that the isotropic 
magneto-electric effect or the so-called 6>-term, Cq = 
{Oe^ /27ih)E ' B with 6> = tt, is suggested to describe the 
the time reversal invarint (TRI) topological band insula- 
tor (TBI) in three spatial dimensions^ (3d). Usually the 
signature of the 3d TBI is the existence of an odd number 
of Dirac cones on the boundary surfaces, when the time 
reversal symmetry (TRS) is preserved^^ When the surface 
states are gapped out by breaking the TRS locally on the 
boundaries, a half integer Quantum Hall effect will take 
place, and give rise to a quantized bulk magneto-electric 
response.^ It is later shown that this isotropic response 
is only a part of the more general anisotropic orbital ME 
tensor defined in the bulk^: 



(1) 



we define ae = |Tra^j and a^d is therefore traceless. 
aij describes either the orbital magneto-polarizability 



(OMP) or the orbital electric susceptibility (OES): 
aij = dPi/dBj = dMj/dEi; 



(2) 



OMP and OES are equal via a Maxwell relation. 

One peculiarity of the ME tensor is that ae is only 
determined up to integer multiples of e'^/h by the bulk 
band structure*^*^ The specific value of ae depends on the 
details at the boundary. From the polarization perspec- 
tive, if we attach an integer quantum Hall (IQH) layer 
with filling v = ±.1 respectively on the top and bottom 
surface of a cylinder, the orbital magneto-polarizability 
(OMP), i.e., dP/dB^ along the axis of the cylinder will 
change by e^//i, due to the density locking to the mag- 
netic field of the top and the bottom IQH layer. From 
the magnetization perspective, if we attach an IQH layer 
with filling v = 1 on the side surfaces of the cylinder, the 
orbital electric susceptibility (OES), dM/dE", will also 
change by e^//i, due to the extra Hall current flowing on 



the surface in response to the electric field. Either way, 
the ME effect is only determined up to an integer multi- 
ple of /h. In the bulk, this ambiguity corresponds well 
to the fact that 6> as a coefficient in front of {E • B) is an 
angle only defined up to integer multiple of 27r, because 
J E • Bd^xdt is quantized, a is odd under TRS, but this 
ambiguity makes it possible for ae not to vanish with 
TRS preserved, as ae = ±e^/2/i is differed from each 
other by /h. ae = ^ and ae = /2h then describe two 
different insulating states of matter under TRS. 

However, the above properties raise some questions. 
From the polarization perspective this ambiguity from 
the bulk is not so surprising, since the zero field po- 
larization is already ambiguous with periodic boundary 
conditions j^i^ From the magnetization perspective, how- 
ever, this ambiguity is a bit more puzzling, because one 
commonly regards magnetization as a bulk property^^i^ 
In particular, with periodic boundary conditions there 
seems to be no reason to expect any ambiguity in the 
magnetization, whereas the ambiguity in the polarization 
is easily understood. 

Before answering this rather specific question, we note 
that there is a even more general one: to what extent 
are thermodynamic quantities such as polarization, mag- 
netization, and ME response determined by the bulk? 
Unlike the conventional thermodynamic quantities which 
are entirely independent of the boundary, we have al- 
ready seen that boundary sometimes plays a role. How 
can we tell when will a thermodynamic variable depend 
on the boundary and when will it not? 

In the following, we will discuss case by case from 
the ground state polarization, orbital magnetization, to 
the magneto-electric tensor. We will argue through 
Gedanken experiments that some of them depends on 
the boundary while others don't. We will verify our ar- 
gument with numerical simulations. By matching the 
observations with our previous calculation done with pe- 
riodic boundary conditions, we can then directly tell from 
the calculation with periodic boundary conditions how 
different thermodynamic quantities depend on bound- 



aries. 
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II. GROUND STATE POLARIZATION 

The ground state polarization is given by the fohowing 
formula with periodic boundary conditions)^ 




In one spatial dimension (Id), the polarization is de- 
fined modulo e with periodic boundary conditions: P = 
Pq + ne, with n an integer. This corresponds to the ob- 
servation that with periodic boundary conditions, we can 
move every electron to the next unit cell and return to 
the original state, while the two states should by defi- 
nition have polarization differed by e. With two ends, 
the polarization will take one specific value, depending 
on the number of charges we put at the two ends. 

However, if there are zero modes at the two ends, the 
polarization is then ambiguous, as theoretically we can 
consider superposition of states of different occupancy of 
the zero modes. The bulk value of the polarization thus 
depends entirely on the boundary. 

In 3d it is a bit more interesting. For simplicity let 
us assume the system sits on a cubic lattice of size a. 
Now the bulk formula has an ambiguity of e/a^, which 
also corresponds well to the fact that we can move ev- 
ery electron to the next unit cell and return to the same 
state. However, with boundary surfaces the situation 
becomes quite different. Consider a capacitor setup. We 
are allowed to put any number of charges on each of the 
opposing surfaces, resulting in a change of the polariza- 
tion in units of e/A (A is the total surface area). In the 
thermodynamic limit, we can put any finite density of 
charges on the surface, and the polarization in the bulk 
can take any value. Our bulk formula is thus no longer 
valid. To accommodate the charge on the surface, how- 
ever, the system needs to either be in a metallic state near 
the boundary, or to break the lattice translation symme- 
try in the two in-plane directions. If neither condition 
is satisfied, then we can only add an integer number of 
electrons per unit cell, and the bulk formula is recovered, 
with the remaining ambiguity determined by the surface. 

How can the bulk formula become invalid? We note 
that the ground state polarization can be understood as 
a Berry's phase when one adiabatically turns on the elec- 
tric field. Firstly, in order for the Berry's phase to make 
any sense, the system has to be gapped. This is the rea- 
son why a metallic surface can render the bulk formula 
invalid. Secondly, if we break the lattice translation sym- 
metry in the two directions perpendicular to the electric 
field, we can no longer integrate over the momentum in 
those directions but should instead sum over a large num- 
ber of sub-bands labelled by the remaining momentum 
along the direction of the electric field. The polariza- 
tion will have an ambiguity of e/A in this case. This 
is different from the conventional thermodynamic quan- 
tity, which will require a symmetry breaking in the bulk 
to change its value. The Berry's phase is thus a rather 



fragile thermodynamic quantity. 



III. GROUND STATE ORBITAL 
MAGNETIZATION 

It is not immediately obvious that the orbital magne- 
tization is independent of the boundary. In the bulk the 
operator M oc (r x v) is ill-defined with periodic bound- 
ary conditions, and seems to be growing as one goes near 
the boundary. Indeed, when one numerically compute 
(M) summing over the local orbitals, there is a finite 
contribution from the boundary orbitals, which renders 
the total orbital magnetization different from the naive 
bulk valueJ^ Nevertheless, it has been shown^^ that the 
boundary contribution is in fact independent of the de- 
tails at the boundary via the use of local Wannier func- 
tions, in an insulator with zero Chern number. 

However, in a Chern insulator, a local Wannier func- 
tion can not be foundiiii^, because the Bloch functions 
cannot be periodic and smoothly defined over the Bril- 
louin Zone. To see that even in this case the orbital 
magnetization is still independent of the boundaries, we 
can consider the following setup: 

Suppose we have an insulator with a non-vanishing 
Chern number in two dimensions. Let us imagine putting 
an auxiliary layer of insulator on top, with an opposite 
Chern number, without any interaction with the origi- 
nal one. The new insulator as a whole is then of total 
Chern number zero. We can therefore make a local Wan- 
nier orbital, by a linear combination of orbitals from the 
two layers The argument then goes through for the in- 
sulator as a whole, and the total orbital magnetization 
should be independent of the boundary. Now since there 
is no interaction between the two layers, the total mag- 
netization is just the sum of the magnetization of the 
original insulator and the auxiliary insulator. We now 
consider a particular boundary condition, where the two 
insulators couples to independent boundary terms that 
do not interact with each other as well. Let us only vary 
the boundary terms that couple to the original insulator. 
The total magnetization cannot change, and neither the 
contribution from the auxiliary insulator. We thus have 
to conclude that even for a Chern insulator, the orbital 
magnetization is independent of the boundaries. 

From this abstract point of view, the generalization to 
Chern insulators seems rather trivial. However, the pres- 
ence of gapless chiral edge states may cause one to worry. 
Suppose we can gate the material to supply a constant 
chemical potential, what will happen if we turn up the 
electric potential on the edge? Will the edge current de- 
crease because fewer edge states are occupied, or will it 
stays the same as required for the bulk magnetization not 
to change? 

We do a straightforward numerical simulation to re- 
solve this paradox. The result is shown in Fig. [TJ We 
can see that while shifting the overall chemical potential 
creates circulating currents, altering the electric poten- 
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FIG. 1. We take our Hamiltonian to be 

^ - J2n ci{Tz-iTx)Cn+x+cl,{Tz-iTy)Cn+y-\-mcl,TzCri + h.C., 

where r's are the Pauh matrices. At half fining with m = 1.5, 
the band carries a Chern number Ci — 1. If we set the chem- 
ical potential /x = 0, the ground state has no magnetization. 
We put the Hamiltonian on a 10 x 10 lattice, and take open 
boundary conditions in both directions. The current on the 
vertical links is plotted. We relate the current to the mag- 
netization by — e^^daM, and take the magnetization at 
the middle to represent the bulk magnetization, (a) /i — 0.5. 
As expected, some edge states are occupied and give rise to 
a bulk magnetization, (b) If we set /j — but locally apply 
an electric potential V = —0.5 to the first two rows at the 
boundary, the edge states are again occupied. However, in 
the region next to those layers, a counterpropagating current 
takes place. The bulk magnetization remains zero (barring 
some finite size effect). 

tial locally at the edge does not change the bulk mag- 
netization. If we look closer, while the current right at 
the edge is changed, there is a counter-propagating cur- 
rent near the edge, which keeps the total current local- 
ized near one edge constant. The counter-propagating 
current is just the integer quantum Hall response to the 
electric potential gradient. This bulk quantum Hall cur- 
rent exactly compensates for the current carried by the 
now-unoccupied edge states, and leaves the bulk magne- 
tization insensitive to the change of the potential local 
near the edge. 

A very similar puzzle arises in the Sz conserved 
spin Hall insulator. On the edge there are counter- 
propagating TR-paired edge states. When we apply a 
uniform Zeeman field Hz, there will be a net circulating 
current from the edge states. We can therefore deduce a 
bulk orbital magnetization response to the Zeeman field. 
We call this the orbital- Zeeman susceptibility. However, 
one can locally break the Sz conservation together with 
the TR symmetry near the edge, to gap out the edge 
states. In this case, will there still be a bulk magnetiza- 
tion response to the Zeeman field? 

The numerical result is shown in Fig. [2j Here we can 
see that even though the edge states are gapped out by 
the local perturbations, the total current flowing near the 
edge remains the same. The local perturbation transfers 
the current from the states at the Fermi level, to the 
occupied bands. In the end, while local properties can 
affect the gapless states, the total current near the edge 
in the is unaffected. 
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FIG. 2. We now think of the previous model as from spin 
up electrons and pair it with its time reversal. We ap- 
plied a unifrom Zeeman field 6Hz — 0.2 cl^SzCn- (a) We 
plot the eigenstate energies in ascending order. The edge 
states live inside the gap. (b) By applying a time-reversal 
as well as Sz symmetry breaking term near the boundary 
6H = Xlneedge cliSxCn, wc cau gap out the edge states, (c)-(d) 
We look at the current on the vertical links. While the cur- 
rent distributes slightly differently with or without the sym- 
metry breaking term at the edges, the contributions to the 
bulk magnetization are identical. 



We therefore conclude that the orbital magnetization, 
as well as the orbital- Zeeman susceptibility is indepen- 
dent of the boundary for an insulator. While the cir- 
culating current may be carried by the edge states, the 
total amount is entirely insensitive to the local boundary 
conditions. One can understand this from a calculation 
with periodic boundary conditions: the magnetization is 
calculated as an energy density in a magnetic field. The 
total energy, unlike the Berry's phase, is a truly extensive 
property, so that the boundary contribution is irrelevant 
in the thermodynamic limit. The energy density in the 
bulk is thus entirely independent of the boundaries far 
away enough, whether there are gapless states or not. 



IV. MAGNETO-ELECTRIC EFFECT 

After the discussion of the polarization and the magne- 
tization and seeing that they are thermodynamic quan- 
tities with very different behaviors, it is thus a natural 
question to ask the same question about the ME ten- 
sor; in addition, about how the Maxwell relation can be 
maintained. Before going into details of the boundary de- 
pendence, however, let us first show that the anisotropic 
part asd is independent of the boundaries. 

In terms of electronic Green's functions and with peri- 
odic boundary conditions, we have derived the ME tensor 
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from the OMP perspective, as a Berry's phase in a mag- 
netic field: 

awzwij = -Y^cibcdT^^^{gdag~^gdbg~^gdcg~^gddg~^g)Sij; 

c^sdij = -^^abjT^^igdiQ'^ gdag~^ gdbg~^ g - h.c). (4) 

The traces include the frequency and momentum integral 
divided by factors of (27r); the symbol Tr*^ denotes the 
integral and trace in one extra dimension in momentum 
space, with the original Brillouin zone and a trivial test 
system as the boundary. While the entire ME tensor 
is derived as a Berry's phase, does not depend on 
the Green's function extended to the extra dimension. 
Without considering boundaries directly, we can show 
that is independent of the boundaries, by showing it 
extends smoothly to finite frequency and momentum. 

At finite frequency and momentum, the ME response 
is understood as a term in the effective action which is 
proportional to E'^{q,uj)B^ {—q, —uj). Unlike the uniform 
ME response however, this term can no longer be un- 
derstood as OMP or OES, due to the fact that unlike 
uniform electromagnetic fields, the electric and magnetic 
fields at finite frequency and momentum are related by 
Faraday's law. The term nevertheless affects properties 
of the propagating electromagnetic waves. For our pur- 
poses, it suffices to show that the effective Lagrangian is 
continuous from g = to g — > 0. At any ^ 0, we can 
calculate the effective Lagrangian by the conventional di- 
agrammatic method. Calculated in the Appendix, the 
bubble diagram gives 

SME = -j -^B'{q)E''{-q)asdk( + 0{q). (5) 

Comparing with Eq. (|4]), we see that asd is continuous, 
whereas o^wzw is entirely absent at finite momentum. One 
might worry that we have missed awzw in momentum 
space due to the fact that it is a total derivative in real 
space, which Fourier transforms to zero and cannot be 
seen in momentum space. However, one can evaluate the 
diagram in real space, and it is still absent. Fundamen- 
tally this is due to the fact that the conventional per- 
turbation theory is perturbative in orders of the gauge 
field, which breaks down with uniform field strength. 
Nevertheless, combining the two calculation, we can still 
say that asd is a bulk property and is independent of 
the boundaries, awzw, on the other hand, is similar to 
the polarization: it does depend on the boundary, but 
when there is no boundary, it presents itself as a Berry's 
phase. Note that one benefit of using the Green's func- 
tion is that the separation of the local terms and bound- 
ary terms matches exactly how the expression depends 
on the extra dimension or not. This is not the case if 
we use the density matrices, either to calculate the same 
Berry's phase^^, or to calculate a current response to a 
pumping procedure^. In both calculations the ME ten- 
sor naturally separates into two terms, with the first term 



independent of the energy gap: 

a = acs + <^G; (6) 

acs is isotropic, but ac is not traceless. While ac can 
be uniquely determined by the bulk band structure and 
is independent of the boundaries, its trace is actually not 
measurable in the bulk. 

Let us now focus at the isotropic part a^zw In terms 
of polarization in a magnetic field, the ambiguity is no 
surprise. However, how does the ambiguity of the orbital 
magnetization in a electric field come about? 

One origin of the ambiguity is from the fact that the 
perturbation of a uniform electric field grows with dis- 
tance. It therefore naturally depends on the boundary, 
when there is one. When we consider periodic boundary 
conditions, however, it becomes less clear. 

In order to study the OES with periodic boundary 
conditions, we first have to properly define the magne- 
tization with periodic boundary conditions. Without the 
current at the boundary, one sensible definition of the 
magnetization is from the relation B = H -\- M. That is, 
in the absence of applied current (which generates i^), 
the magnetization simply equals the measured magnetic 
field. Note that with periodic boundary conditions and 
a finite volume, the magnetic field is quantized, because 
the total magnetic flux through the sample is quantized 
in units oih/e. In this case we take the perspective that 
the magnetic field will take the closest quantized value to 
the magnetization while the magnetization itself is still 
continuous. 

In our previous worki^, we have shown that in a mag- 
netic field, the 6 term, which characterizes the isotropic 
part of the OMP, changes the quantization condition of 
the global electric flux. The ground state of the sys- 
tem thus carries an electric flux of -{Oe'^ /27rh)^B + ^e, 
where n is some integer that minimizes the flux. Using 
= D = E-\-P^ the 6 term thus gives an isotropic orbital 
magneto-polarization response = However, this 
result is valid only when (<l>^6>e/27r/i) < 1. In the ther- 
modynamic limit this condition is always violated, and 
instead f§ = 0. 

Similarly, to see whether the same term contributes 
to the OES of the system, we would like to investigate 
whether there is a uniform magnetic field, when we con- 
strain the path integral to have a given average electric 
field in the same direction. However, the electric field and 
the magnetic field behave in intrinsically different ways, 
when we formulated our theory assuming the existence of 
electric charges and the absence of magnetic monopoles: 
the quantization of the electric flux can change in the 
presence of the magnetic field, while the quantization of 
the magnetic flux is fixed at ih/e). When we apply an 
electric flux, we can always imagine that the system is a 
coherent state composed of states with integer electric 
fluxes. The background magnetic field therefore does 
not have to be different from zero. Therefore, even at 
finite size, the term does not give rise to the OES. The 
Maxwell relation between the the isotropic OMP and the 
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OES are thus violated. They are only equal in the ther- 
modynamic limit, where the term gives no contribution 
for both quantities. In other words, the isotropic OES 
is better thought of as a bulk-induced surface response, 
which vanishes when there is no boundary surfaces. 

Now let us consider geometry with boundaries in some 
detail. From the result of Ref.^, we know that with open 
boundary conditions in all directions, the OES has an 
ambiguity only determined by specific surface boundary 
conditions. We have also seen in the introduction that in 
a cylinder geometry, the ambiguity of the OES can come 
from the quantized Hall current on the side surfaces. 




0.25 



the magnetic field will converge to the expected value, in 
contrast to the situations with periodic boundary condi- 
tions where it stays at zero. We have numerically con- 
firmed this result by calculating the magnetization in the 
electric field, using the momentum space formula for the 
magnetization, derived in Ref^ii^, as shown in Fig. [3l 

Before summing up, let us consider how gapless surface 
states can alter the ME response. Evidently, if we attach 
a fractional quantum Hall state on the side of the cylin- 
der, the OES is going to change by a fraction of e'^/hi^ 
In general the fraction is quite arbitrary, so in this case 
the bulk value of the isotropic OES is not valid. This 
corresponds to the fact that the fractional quantum Hall 
state has ground state degeneracy. In general, we will 
therefore expect any gapless surface state will destroy 
the bulk description of the isotropic ME response. 

To sum up. The anisotropic part of the ME tensor 
is independent of the boundaries. The isotropic part 
o^wzw depends partially on the boundary. While asd is a 
truly local quantity, awzw only lives at g = 0. Corrobo- 
rating with the fact that both isotropic OES and OMP 
responses vanish with periodic boundary conditions in 
the thermodynamic limit, it is better to think of 

C^wzw ^S 

a quantized surface effect induced by the bulk. 



FIG. 3. Here we plot the calculated OES versus the number of 
layers in the 2;-direction, with the model described by Eq. (73) 
in Ref.- with 6' = 0, m = c= l. We take O.bn at the 
top and the bottom layer to gap out the edge states. (If we 
take = =b0.57r on the two surfaces respectively, the whole 
system will be a Chern insulator and can no long be kept at 
charge neutrality without closing the gap in a magnetic field.) 
We put on an electric field such that the potential difference 
between the top and the bottom layer is 0.2. The boxes show 
the calculated values. The solid curve is a fit by assuming a 
fixed width w of the surface charges when there is a magnetic 
field, such that ^ oc (1-f )• The fit gives gives ^ = 0.50^ 
in the thermodynamic limit and w = 2.54. The OES changes 
sign as expected, when we change to = — O.Stt instead on 
the boundary. 

What if there are no side surfaces? Suppose we take pe- 
riodic boundary conditions only in two directions to get 
rid of the side surfaces. Does the OES still have the same 
ambiguity? One naively would expect the situation to be 
similar to the case with periodic boundary conditions, 
due to the absence of the possible circulating Hall cur- 
rents. However, a more careful argument shows it is not 
the case. In fact, the system will spontaneously generate 
a magnetic field, which will then generate surface charge 
density a = ±{i'-\-6/27r)e'^B /h via the OMP response, to 
lower the electric energy. Minimizing the total energy as 
a function of B, we then get B = M = {u -^0 /27i)e^E /h. 
While at finite size the total magnetic flux is quantized 
in units oih/e'm this setup, in the thermodynamic limit. 



V. CONCLUSION 

We have thus gone through polarization, magnetiza- 
tion, and magneto-electric responses and see their de- 
pendence on the boundary. A bulk calculation done with 
periodic boundary conditions contains enough informa- 
tion to predict how the quantity in question can depend 
on the boundary, however. In particular, using our for- 
malism described in Ref.^^, any quantity that does not 
involve an extension of the Green's function to one ex- 
tra dimension is independent of the boundary. On the 
other hand, quantities that requires an extension to ex- 
tra dimension will depend on the boundary. The bulk 
can determine its value up to some quantized amount, 
only when (i) there are no gapless surface states, (ii) sur- 
faces break no symmetry that is required to determine 
the bulk value with periodic boundary conditions, and 
(iii) the system is kept at charge neutrality. If any of 
the conditions are violated, the surface contribution will 
dominate and render the results obtained with periodic 
boundary conditions invalid. 
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APPENDIX: ME EFFECT AT FINITE MOMENTUM 



In the main text we have argued heuristically that the trace of the ME tensor comes entirely from the surface, and 
therefore does not contribute at finite momentum. The locally measurable ME tensor is therefore traceless in the 
q ^ limit. We can directly calculate the ME tensor at finite q: Fourier transforming and expanding the hopping 
Hamitonian up to second order of A^, we have 



k^q k,q,q' 

Hk = ex.p{ikdi) is a matrix. Integrating out the electrons, the effective action at quadratic order of reads 

similarly, the trace includes the integral of energy and momentum divided by 27r. The first term in the trace is from 
the second term in Eq. (O, usually called the paramagnetic current, and does not have q dependence. To compare 
with Eq. (|4]), we Taylor-expand the second term to second order in q to get the behavior in the q ^ limit: (From 
here on, we drop the dependence on k to avoid cluttering.) 

(9) 

To further simpHfy the expression, let us now take the Coulomb gauge. In the Coulomb gauge, we have to take 
either A or <j to be in the time direction to have the expression contribute to the ME tensor. Since didojg~^ = 0, we 
can integrate-by-part the time derivative, using d^jQ = —g^^ and rename the indices j, /c, now running through only 
the spatial directions, we get 

Sme^ I -^'-i^q'A^q)A''i-q)Tv(^gd,g-'d,gdkg-'gy (10) 

Now we need to massage the expression a little bit. Let us use (ijk) as a short-hand notation of the expression 
TrI gdig~^djgdkg~^g )• Integrating by partsi^, we have the following relation: 



{ijk) + {jki) + {kij) = 0. (11) 

We therefore have 

(jik) = ^ {2{jik) - (ikj) - (kji)). (12) 
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In the trace in Eq. ([TOj), only the part symmetric under the exchange of the index j and k would contribute, as we 
can change variables from q to — effectively exchanging (q) and A^{—q). Therefore, in the expression above, we 
can exchange j and k freely. We therefore have 

Sme^ I ^lujq^A^{q)A''{-q){{jik) + {kij)-{ijk)-{kji)) 

/d^ i 
j^^Loq^A^q)A''{-q)€,jeeabe{{kab) + ihak)). 

= j ^^'-B\q)E\-q)eM{{kah) + {bak)) 

^ - / -^^B\q)E'^{-q)au{q ^ 0). (13) 

aki{q ^ 0) = —^eau{{kcib) + ihak)) is traceless, as the two terms cancel each other with antisymmetrization. 
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